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Multipartite Bell-type inequalities are derived for general systems. They involve up to eight 
observables with arbitrary spectra on each site. These inequalities are closely related to the algebras 
of quaternions and octonions. 
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In their famous paper Einstein, Podolsky and Rosen 
(EPR) suggested a Gedankenexperiment [l[ . As they be- 
lieved, it would prove the incompleteness of quantum me- 
chanics. An interesting analysis of this problem was given 
by Bohr Q , who did not agree with EPR. Almost thirty 
years later, in 1964, in his remarkable paper Bell pro- 
posed a quantitative test which should resolve the prob- 
lem of completness of quantum mechanics 0] ■ He showed 
that the assumptions of the EPR arguments lead to some 
restrictions on multipartite correlations, which are now 
referred to as Bell inequalities. If Bohr's arguments are 
correct, then these inequalities can be violated. In the 
1980's the first experimental violations of these inequali- 
ties were demonstrated 0, @l ; an -d thus the arguments of 
Bohr were verified. 

There exist many different Bell- type inequalities. Most 
of them deal with dichotomic observables or with gener- 
alizations to observables with more general discrete spec- 
tra, for example see |6j . On the other hand, in the original 
EPR paper the situation of observables with a continu- 
ous spectrum was considered. For the case of observables 
with a general spectrum, to this end the theory of Bell- 
type inequalities is much less developed. 

The first multipartite Bell-type inequality, valid for ar- 
bitrary observables, has been obtained very recently 0]. 
In the simplest, bipartite, case with two observables Aj, 
Bj on each site (j = 1, 2) it reads as 



K(ii +iB 1 )(A 2 + iB 2 ))\ 2 = (A 1 A 2 -B 1 B 2 f 
+ (BiA 2 + A 1 B 2 ) 2 < ((If + Bl){A\ + B 2 2 )). 



(1) 



The original proof is based on ignoring local commuta- 
tors in the quantum mechanical analogue of a classical 
inequality. In some sense such a procedure is ambiguous, 
since only those commutators appearing explicitly are ig- 
nored. Furthermore, no direct proof was given that the 
resulting inequality is fulfilled for any separable quantum 
state - only in this case it should be called Bell-type. 

In this letter we give a strict and unambiguous proof 
of general Bell-type inequalities, including their relation 
to the separability problem. It allows us to consider up 
to eight arbitrary observables at each site, for a com- 
plex quantum system. The multipartite extension of the 
inequalities is based on the algebra of quarternions and 
octonions. The recently proposed inequalities will ap- 
pear as special cases of our approach. 



Let us start to prove that each bipartite separable state 
satisfies the inequality |T]). This approach also allows 
us to obtain Bell-type inequalities with four and eight 
observables on each site. First, we prove a more general 
statement: if F\, . . . , F% and G are Hermitian operators 
such that some (in general, multipartite) states Qj, j — 
0, 1, . . ., satisfy the inequality 



(F 1 ) 2 + ... + (A> 2 < (G), 



(2) 



then arbitrary mixtures (i.e. convex combinations) of 
these states also satisfy this inequality. The proof is 
based on the Cauchy-Schwarz inequality 



IE' 



(3) 



which is valid for all complex numbers Xj and yj. The 
case of yj = pj being a probability distribution, i.e. pj ^ 
and Ylj=^>Pj = 1, is of special interest. Then one gets 
from the inequality (|3j) 



IE^" 5 

j=0 



3=0 



(4) 



Note that the sum on the left hand side of this inequality 
is the mathematical expectation of the random variable 
X that attains the value Xj with probability pj , and the 
right hand side is the mathematical expectation of the 
square of this random variable. The inequality ([4]) then 
states that variance of X, ax, is nonnegative. 

Let us now take a convex combination g = X) jiTo Pj &j 
and estimate the left hand side of the inequality @ . Ac- 
cording to the inequality (fJl) we have 



(Ft) 2 



(T,Pj{Fl)j 



E 

j=0 



Pj(Fi) 



(5) 



where (Fi)j is the average value on the state Qj. Now we 
can estimate the left hand side of ([2]) as follows: 



k +oo k +oo 

^<f ; ) 2 < $>E<% 2 ^ = ( 6 ) 

j=0 1=1 3=0 



1=1 



So, we have obtained the desired result. 
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Now we can easily prove the inequality (fT]) for separa- 
ble states. Here we have k = 2 and F\ — A\A 2 — B1B2, 
F 2 = B X A 2 - A 1 B 2 and G = (if + B 2 )(A 2 2 + B\). The 
have to show that the inequality |T]) is valid for all fac- 
torizable states. Then we get 

(A ± A 2 - B\B 2 ) 2 + (B X A 2 + A X B 2 ) 2 = 
((A,) 2 + (B,) 2 )^) 2 + (B 2 ) 2 ) < (7) 
((A 2 + B 2 )(A 2 +B 2 )). 

The equality in this chain is valid since only squares re- 
main in the sum and mixed terms cancel each other. The 
last step, the inequality, is valid since it just expresses the 
fact the variance of an observable is nonnegative. We con- 
clude that each convex combination of factorizable state, 
i.e. each separable state satisfies the inequality fT]). 

We see that the key point in our proof of the inequality 
(fT]) is the estimation ([7]), which can be divided into two 
steps. The first step, the equality, can be expressed as 
the following square identity: 

(01 a 2 - M2) 2 + (0162 + ha 2 ) 2 = (af + b\){a 2 2 + b\). (8) 

The second step is valid due to the nonnegativity of the 
variance of observables. Having an identity of the form 

(a 2 + b\ + .. .){al + b 2 + ...)= x 2 + y 2 + ... , (9) 

where all the sums contain the same number n of terms 
and x, y, ... are bilinear functions of ai, bi, . . . , I = 1, 2, 



we can immediately write a Bell-type inequality 

(X) 2 + (Y) 2 + . . . < ((A 2 +B 2 + . . .)(A 2 2 +B 2 2 +. . .)), (10) 

where X, Y, . . . are the Hermitian operators obtained by 
replacing a; , bi, ... by arbitrary Hermitian operators A[ , 
Bi, . . . , I — 1, 2 in the bilinear forms x, y, . . . respectively. 
Let us again formulate the reasons why this inequality is 
valid for all separable states. Firstly, it is of the form 
i.e. if it is valid for some states, it is also valid for their 
mixtures. Secondly, it is valid for all factorizable states 
due to the identity ((9J and nonnegativity of the variances 
of observables. It follows that it is valid for all mixtures 
of any factorizable states, i.e. for all separable states. 

What square identities exist? The case of n = 2 was 
considered above. There are also square identities for 
n = 4 and n — 8. They read as: Euler four square 
identity 

(01O2 - bib 2 - C1C2 - did 2 ) 2 + 
(61 02 + aib 2 - dic 2 + cid 2 ) 2 + 
(cid2 + dib 2 + aic 2 - bid 2 ) 2 + (11) 
(dia 2 - C1&2 + b\c 2 + aid 2 ) 2 = 
(4 + b\ + c\ + d 2 )(a 2 + b 2 2 +c 2 2 + d 2 ), 

and Degen eight-square identity 



(0102 - bi&2 - C1C2 — diGb — eie2 - /1/2 - 3132 - /u/12) + (6102 + ai&2 + diC2 — cicb + fiez - ei/2 - higz + 31/12) + 
(cia 2 - dib 2 + aic 2 + bid 2 + 3ie 2 + /11/2 - ei3 2 - /i^2) 2 + (dia 2 + C162 - bic 2 + aid 2 + feie 2 - 31/2 + /132 - ei/i 2 ) 2 + 
(eia 2 - /162 - 31 c 2 - hidi + aie 2 + 61/2 + C132 + dih 2 ) 2 + (/1O2 + e±b 2 - h\c 2 + g±d 2 - b±e 2 + a±f 2 - dig 2 + c 1 h 2 ) 2 + (12) 
(31O2 + hib 2 + eic 2 - j\d 2 - ae 2 + dif 2 + a\g 2 - bih 2 ) 2 + (hiaz - g\b 2 + /1C2 + e\d 2 - d\e 2 - C1/2 + fri3 2 + a±h 2 ) 2 = 
[a\ + b\ + c 2 + d\ + + f 2 + g\ + h\)(a 2 2 + b\ + c 2 + d 2 2 + e\ + /| + 32 + hi). 

I 



The famous Hurwitz theorem states that there are no 
other identities of such a form Q. 

The inequality corresponding to the identity (fTTj) is 

(ilia - B X B 2 - C X C 2 - D 1 D 2 ) 2 + 
(Bii 2 + A X B 2 - D1C2 + C 1 D 2 ) 2 + 
{C 1 A 2 + D 1 B2+A 1 C2-B 1 D 2 ) 2 + (13) 
(DxA 2 - C X B 2 + B X C 2 + A 1 D 2 ) 2 < 
((if + B 2 + C 2 + D 2 )(A 2 + B 2 + C 2 + D 2 )) 

This is a bipartite Bell-type inequality with four observ- 
ables on each site. To extend it to the general multi- 
partite case it is natural to use the algebra of quater- 



nions. Remember that this algebra has dimension 4 
over the reals, so each quaternion q can be written as 
q = x + iy + ju + kv in an unique way, where x, y, u and 
v are reals. The multiplication rules for the imaginary 
units i, j and k are i 2 = j 2 = k 2 = — 1, ij = —ji = k, 
jk = —kj — i, ki — ~ik = j . The conjugation q* of the 
quaternion q is defined as q* = x — iy — ju — kv. The 
norm of q is defined in the standard way as \q\ = \Jq*q. 
The identity pT|) simply express the fact that the norm 
is multiplicative: \q'q"\ = \q'\\q"\ for arbitrary quater- 
nions q' and q" . The norm also satisfies the triangle 
inequality \q% +q 2 \ ^ |<7i| + \qz\- The inequality (fT]) is for- 
mulated using the operators of the form / = A + iB 1 
which is a general form of a non-Hermitian operator. 
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We extend the class of operators acting on the state 
space of the system to quaternionic operators of the form 
q = A + iB+jC+kb, where A, B, C and D are ordinary 
Hermitian operators. Since the algebra of quaternions is 
noncommutative, care must be taken when defining the 
product of quaternions with operators. We define this 
product such that, if q m = A rn + iB rn + jC m + kD m 
are quaternionic operators acting on different degrees of 
freedom m — 1, . . . , n, then 

{q x . . . q n ) = (ft) . . . <g„) (14) 

for each completely factorizable state. Let us calculate 
the average value of the product jf, where the operator 
/ has been defined above 

07") = .l(f) = J((A) + i(B)) = j(A) ij{B) 

(15) 

= (A — iB)j = (Pj). 

Here we used the fact that the numbers (A) and (B) are 
real and due to this they commute with j. The same 
is valid with respect to the other imaginary unit k. We 
see that the natural way to define the product of the 
quaternion q with the operator / is qf = f(x + iy) + 
P (ju + kv). In particular, if F is a Hermitian operator, 
then we have qF = Fq, so the quaternionic operators q 
defined above behave as if the operators A, B, C and 



D were ordinary real numbers. This guarantees that the 
equality (fT4|) is fulfilled for all factorizable states. 

The same idea we used to prove the inequality ([2]) al- 
lows us to prove the following statement: if q m are quater- 
nionic operators and F m are Hermitian operators, acting 
on different degrees of freedom, such that | (q m ) | 2 < (F m ), 
m = 1, . . . ,n, then each completely separable state sat- 
isfies the inequality 

\(q X ...q m )\ 2 ^(F x ...F n ). (16) 

Here we need the factorization property (|14p , multiplica- 
tivity of the norm and the triangle inequality. Since we 
can estimate |(qm)| 2 as 

\(q m )\ 2 = (A m ) 2 + (B m f + (C m ) 2 + (An)' 

< (A 2 m + B 2 m + C 2 n + D 2 m ), 

we can take F m — A 2 m + B^ + + b 2 m . Upon taking 
the product of quaternionic operators q\ , . . . , q n , the left 
hand side of the inequality ()16|) will be a sum of four 
squares of average values of some observables. Then the 
inequality (|16p is a multipartite Bell-type inequality with 
four observables on each site. In the case of n = 2 it is 
the inequality (fT3|) . 

The inequality corresponding to the identity (fT2|) is 



(A1A2 — B1B2 — C1C2 ~ D1D2 — E1E2 — F1F2 — G1G2 — H\H2) 2 + 

(B1A2 + A1B2 + D1C2 — C\F>2 + F\Ej2 — F:\F2 — H1G2 + GiJ?2) 2 + 

(C 1 A 2 - b x B 2 + A&2 + B X D 2 + G X E 2 + H X F 2 - E X G 2 - F 1 H 2 ) 2 + 
(D1A2 + C1B2 — B1C2 + A\F>2 + H\Ej2 — G\F 2 + F\G 2 — EiH 2 ) 2 + 

{E 1 A 2 - F X B 2 - G 1 C 2 - HxD 2 + A X E 2 + B^ + C X G 2 + D 1 H 2 ) 2 + (18) 
(F 1 A 2 + E X B 2 - HiC 2 + G x b 2 - B X E 2 + A 1 F 2 - + C 1 H 2 ) 2 + 

(GxA 2 + HiB 2 + E X C 2 - F x b2 - C X E 2 + b x F2 + A X G 2 - B 1 H 2 ) 2 + 
(HiA 2 - G X B 2 + F X C 2 + E x D 2 - D X E 2 - C X F 2 + B X G 2 + A X H 2 ) 2 < 
({A 2 + B 2 + C\ + D 2 + E 2 + F 2 + Gj + H 2 )(A 2 + B 2 + C\ + D 2 + E 2 + F 2 + G\ + H 2 )). 

I 



It is a bipartite Bell-type inequality with eight observ- 
ables on each site. To get a general multipartite inequal- 
ity we need to use the algebra of octonions. It is an 
8-dimcnsional algebra over the reals, so each octonion o 
can be written as 

O = Xo + X 1 i 1 +X2l2 + X3l3 + X4 : i4 + X 5 i 5 +X 6 i 6 + X7i7 (19) 

in an unique way, where xi, I — 0, . . . , 7 are reals and 
il, I = 1, . . . , 7 are imaginary units, whose multiplication 
rules are given by Table |U The conjugation o* of the 



octonion (fl9|) is defined as 

o* = x Q -xiii-X2i2-X3i3-x i i4 : -x 5 i 5 -x 6 i 6 -x 7 i 7 . (20) 
The norm |o| is also defined in the standard way as 



|o| = \fo*o = y xl + x\ + x\ + x"l + x\ + x\ + x\ + x 2 . 

(21) 

The identity (fT2")l expresses the fact that this norm is 
multiplicative: \o'o"\ = |o'||o"|. This norm also satisfies 
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reads as 



TABLE I: The multiplication table of imaginary units of oc- 
tonions. 



the triangle inequality \o\ + 02I 
define octonionic operators as 



< |oi| + |o 2 |. We can 



o = A + iiB + i 2 C + i a D + i 4 E + i 5 F + i 6 G + i 7 H, (22) 

where A, . . . , H are Hermitian operators. Here we iden- 
tify the first imaginary unit i\ with the standard complex 
unit i. For the product of the other imaginary unit with 
ordinary operators / we have the relation: iif — f'ii, 
I = 2, ... ,7. But now the product is not associative, so 
whenever we deal with a product of more than two terms 
we have to explicitly group the terms. 

Now we can generalize the inequality §TE§ as follows: 
if 6 m are octonionic operators and F m are Hermitian op- 
erators acting on different degrees of freedom such that 
|(6 m )| 2 ^ {Fm}, m = 1, ■ . ■ ,n, then each completely sep- 
arable state satisfies the inequality 

\(o 1 ...d n )\ 2 sC (F 1 ...F n ), 



(23) 



for all C n possible groupings of the terms on the left 
hand side, where C n is the nth Catalan number defined as 
C n = (2n — 2)\/{n\(n — 1)!) (we did not show the brackets 
explicitly). Here we can take F m = A 2 m + B 2 n + C m + 
D 2 n + Ef n + F? n + G 2 m + H 2 n . Upon taking the product 
of 01, . . . , d n , the left hand side of the inequality (|23|) 
will be a sum of eight squares of average values of some 
observables. Then the inequality ([23| is the multipartite 
Bell-type inequality with eight observables on each site. 
In the case of n — 2 it reduces to the inequality (fT5|) . 

It is noteworthy that the inequalities |T]), (fTB")) and (fT5]) 
form a hierarchy — the inequality ([T]) is a special case of 
(|13JL which is in turn a special case of (fl8|) . As shown 
in [21], the inequality ([T]) can be violated. Thus it is clear 
that our inequalities (fT3|) and (fT8|) can be violated as well. 
Any violation of these inequalities is a clear signature of 
entanglement. 

The inequalities ([1]) , (fl3|) and (fT8|) can be also obtained 
in another way. The integral form of the inequality (j4]) 



p(X)X(X)d\ 



p(\)\X{\)\ 2 d\ (24) 



where p(A) is a probability distribution on a measurable 
set A, and X{\) is a real- or complex-, quaternion- or 
octonion-valued function of A e A. The set A can be 
thought of as a set of hidden variables, which completely 
specify the state under study. The inequality sim- 



ply states that |(A")| 2 ^ (|A"| 2 ). Let us take instead of 
X the operator X, which is a product of ordinary non- 
Hermitian operators f m , quaternionic operators q m or oc- 
tonionic operators o m defined in pip , acting on different 
degrees of freedom m = 1, . . . , n, and find the quantum 
mechanical analogue of the quantity |A^| 2 . For example, 
in the case of ordinary non-Hermitian operators, we have 
flif™, = -A^ + B 2 ^ + i [A m , B m ] . In a local hidden variable 
theory all commutators must be zero, so \X\ 2 must be re- 
placed by the product Y\^ n= i{A 2 n + B 2 n ). Analogously, in 
the case of quaternionic operators |A"| 2 must be replaced 
by ]lm=l(^m + Bm + ^ra + ^ro) and in tlie case of oc " 
tonionic case by Ilm=i(^m + • • • + H 2 n ). Then we again 
arrive to the inequalities ([1]), (fL3|) and (fLS]) respectively. 
This approach is simpler then the one we started with, 
but it is ambiguous and does not give strictly relate the 
obtained inequalities to separability. 

In conclusion, we have obtained Bell-type inequalities 
for observables with a general spectrum. They apply to 
measurements of up to eight observables, for arbitrary 
systems. The multipartite forms of the inequalities are 
related to the algebras of quaternions and octonions. 
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